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Abstract 

Equivalent crystal theory (ECT) is applied to the study of multilayer relaxations and 
surface energies of high-index faces of Fe, Al, Ni and Cu. Changes in interplanar 
spacing as well as registry of planes close to the surface and the ensuing surface energies 
changes are discussed in reference to available experimental data and other theoretical 
calculations. Since ECT is a semiempirical method, we investigate the dependence of 
the results on the variation of the input used. 
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1. Introduction 


In the last ten years, there has been a large number of experimental and theoretical 
studies on the subject of surface structure of high-index faces of metals [1-13]. Several 
low-energy electron diffraction (LEED) and high-energy ion scattering (HEIS) experiments 
provided a wealth of information, supplemented by numerous theoretical models, on the re- 
laxation and energetics of unreconstructed metallic surfaces which, in the case of high-index 
faces, includes the relaxation of interlayer registries as well as interlayer spacings, where a 
change in the surface-parallel components of the interlayer vectors can occur without di- 
minishing the symmetry of the surface space-group [2]. However, the number of studies on 
high-index faces has been limited to just a few systems (Fe(210),(211),(310) [3]; Al(210) 
[4], (311) [5], (331) [2] and the (311) face of Cu [6] and Ni [7]) thus providing limited infor- 
mation for extracting definite relaxation patterns and global behavior of such properties. 

That is not the case for low-index faces: as summarized in ref. 1, several experimental 
techniques and almost all of the theoretical methods available have been devoted to the 
study of such systems. Although there is yet no simple model that fully accounts for 
the observed patterns in the high-symmetry cases, the available data, both experimental 
and theoretical, gives a good but necessarily incomplete description of general trends and 

patterns. 

We start by discussing the (311) face of fee metals, the fourth highest density surface, 
which has a relatively wide open structure. High index planes have rough contours, for which 
a smoothing of the electron density is likely to occur. Such density smoothing provides 
a driving force for interlayer registry shifts. This same argument helps to understand 
the relatively large contractions of the first interlayer spacing, Ad 12 , as compared to the 
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corresponding values found for low-index faces, in agreement with the observation that Ad u 
increases with surface roughness [3]. 

In a previous paper [12] we focused our attention- on one case, A1 (210), which we 
used both as a testing ground for the application of the theoretical framework provided by 
equivalent crystal theory (ECT) [13], as well as an interesting example on which we based 
a new concept generalizing the idea of roughness of a surface. Because of the excellent 
agreement with experiment found for that system, in this work we conclude our survey of 
multilayer relaxation studies by examining other fee (210) surfaces with equivalent crystal 
theory, for those cases for which experimental data exist. We complete our study with a 
discussion of the structure of fee (331) and bcc (210) and (310) surfaces. 

All the reliable experimental data comes from analysis of LEED experiments. Even in 
this case, the values for relaxations for interlayer spacing and registry come from a rather 
complicated procedure. The relaxations are determined from multivariable, least squares, 
fits to intensity-energy curves for various beams, based on predictions from multiple scat- 
tering theory. Error estimates for the fits are based on assuming a quadratic form for 
deviations from the optimum r- value [4] in the minimization search in terms of the input 
parameters (i.e., layer spacing and registry, potential parameters and Debye temperature). 
The error bars attached to the experimental results for the relaxations thus included infor- 
mation from changes in the input parameters of the multiple scattering model as well as 
errors inherent to the experimental technique used. In this work we address a similar issue, 
by defining the uncertainties in the optimum relaxation values obtained from an energy 
search, by including a margin of error in our theoretical predictions that, in the particular 
case of semiempirical methods such as ECT, can be attributed to fluctuations in the input 
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data used (generally experimentally determined). 

The paper is organized as follows: in section 2 we briefly discuss equivalent crystal 
theory and provide the essential working equations. Section 3 focuses on the case of per- 
pendicular relaxations of low-index faces in order to illustrate the need for the introduction 
of theoretical ‘error bars’, thus providing a better framework for comparison between ex- 
perimental and theoretical results. In section 4 we apply ECT to several fee high-index 
metallic surfaces and compare with experimental values when available. Conclusions are 
drawn in section 5. 

2. Equivalent Crystal Theory 

Equivalent crystal theory [13) is based on an exact relationship between the total energy 
and atomic locations and applies to surfaces and defects in both simple and transition 
metals as well as in covalent solids. Lattice defects and surface energies are determined 
via perturbation theory on a fictitious, equivalent single crystal whose lattice constant is 
chosen to minimize the perturbation. The energy of the equivalent crystal, as a function of 

its lattice constant is given by a universal binding energy relation [14]. 

\ 

Let e be the total energy to form the defect or surface, then £ = £. £>' where £,• is the 
contribution from an atom i close to the defect or surface. ECT is based on the concept 
that there exists, for each atom i, a certain perfect, equivalent crystal with its lattice 
parameter fixed at a value so that the energy of atom i in the equivalent crystal is also £,. 
This equivalent crystal differs from the actual ground-state crystal only in that its lattice 
constant may be different from the ground-state value. We compute £,• via perturbation 
theory, where the perturbation arises from the difference in the ion core electronic potentials 
of the actual defect solid and those of the effective bulk single crystal. 
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For the sake of simplicity, the formal perturbation series is approximated by simple, 
analytic forms which contain a few parameters, which can be calculated from experimental 
results or first-principles calculations. Our simplified perturbation series for £; is of the form 

£| - = A E If* K(t)] + KO'J)] + E k )] + J2 F ' 1 (1) 

{ j i,k ™ ) 

where F * [a*] = 1 - (1 + a*)e _a ’ and A E is the cohesive energy. In this expression, we 

distinguish four different contributions to the energy of atom i and thus, the existence of 

four different equivalent crystals which have to be determined for each atom i. 

The first term, F*[a^(i)], contributes when average neighbor distances are altered 
via defect or surface formation. It can be thought of as representing local atom density 
changes. In most cases, this ‘volume’ term is the leading contribution to a and in the case 
of isotropic volume deformations, it gives £, to the accuracy of the universal energy relation 
[14]. The value of aj(i), the lattice parameter of the first equivalent crystal associated with 
atom i, is chosen so that the perturbation (the difference in potentials between the solid 
containing the defect and its bulk, ground-state equivalent crystal) vanishes. Within the 
framework of ECT, this requirement translates into the following condition from which aj(i) 
is determined: 

NR^e' 0 ^ + MR>e-l a+ ^ R > - £ = 0 ( 2 ) 

def ect 

where the sum over the defect crystal or surface is over all neighbors within second-neighbor 
(NNN) distance. Tj is the actual distance between atom i and a neighbor atom j, N and 
M are the number of nearest-neighbor (NN) and next-nearest-neighbors, respectively, of 
the equivalent crystal (12 and 6 for fee, 8 and 6 for bcc) and p,o and X are parameters 
known for each atomic species, listed in Table 1. S(r : ) is a screening function and R 1 


5 



and R 2 are the NN and NNN distances in the equivalent crystal. The equivalent lattice 
parameter, ai, is thus related to the scaled quantity aj via = (-j 1 - ^Se) /*> w ^ ere 
twse is the equilibrium Wigner-Seitz radius, l is a scaling length and c is the ratio between 

the equilibrium lattice constant and t\yse - 

The higher-order terms are relevant for the case of anisotropic deformations. The 
linear independence attributed to these four terms is consistent with the limit of small 
perturbations which we assume for the formulation of ECT. The second term, F m [a^(i, j)], 
is a two-body term which accounts for the increase in energy when NN bonds are compressed 
below their equilibrium value. This effect is also modeled with an equivalent crystal, whose 
lattice parameter is obtained by solving a perturbation equation given by 

A JR p e~ aRl - NR p 0 e- oH<> + A 2 R P 0 " Ro)e~ m} ~ I ^ ) = 0, (3) 

where /? = 4a for the metals used in this work, and R x is the NN distance of the equivalent 
crystal associated with the deviation of NN bond length R } from Rq, and R 0 is the bulk 
NN distance at whatever pressure the solid is maintained (generally, R 0 is the ground- 
state, zero-pressure value). A 2 is a constant determined for each metal (see Table 1 for 
a list of values of A 2 used in this work). The scaled equivalent lattice parameter is then 

Oj = - WSfi) /*• 

The third term, F* [<£(», j,*)] accounts for the increase in energy that arises when 
bond angles deviate from their equilibrium values of the undistorted single crystal. This 
is a three-body term and the equivalent lattice parameter associated with this effect is 
obtained from the perturbation equation 

NR p e~ aRl - NR p 0 e-° Ro + A 2 Rle~ a{R ^ Rk ~ 2R<,) sin (9 jk - 6) = 0 (4) 
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where A 3 is a constant listed in Table 1 and 0 jk is the angle between the NN distances 
Rj and R k with the atom i at the center. 6 is the equilibrium angle, 70.5 degrees for bcc 
and 90 degrees for fee. This term contributes only when there is a bond-angle anisotropy 
(0 jk ^ 0). The scaled lattice parameter is then a$ = - t W se) H- 

The fourth term, F‘ K(i,p,q)], describes face diagonal anisotropies (see Ref. 13 for a 
detailed description, for each lattice type, of the structural effect associated with this term). 
The perturbation equation reads 

NRje~ aRl - jvsje-* + = 0 (5) 

where d is the face diagonal of the undistorted cube and A 4 is a constant adjusted to 
reproduce the experimental shear elastic constants (Table 1). Finally, a 4 — ^ c f ivse) /^- 
Consider a rigid surface (i.e., no interlayer relaxation): all bond lengths and angles 
retain their bulk equilibrium values, thus F*(a * 2 ) = F'(a|) = F*(a*) = 0. The surface 
energy is therefore obtained by solving for the ‘volume’ term represented by F\a\) only. 
If we consider a rigid displacement of the surface layer towards the bulk, as is the case in 
most metallic surfaces, the higher-order terms become finite: some bonds are compressed, 
contributing to F’(a’ 2 ), the bond angles near the surface are distorted as well as the differ- 
ence between face diagonals in some cases, generating an increase of energy via F*(a$) and 
F*(a’ 4 ). For the case studied in this work, these additional contributions to £, are generally 
small, representing only 1 % to 2 % of the total energy. However, while these anisotropy 
terms are small for metals when there is no reconstruction, they play an important role in 
the energetics of these defects where the differences in energy between the rigid and relaxed 
configurations are also small. 
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3. Multilayer relaxation of pure crystals 

Before proceeding to the calculation of multilayer relaxation in high-index faces, we 
will discuss some features of theoretical calculations of these quantities. Ref. 1 provides 
a reasonably large sample of both experimental and theoretical results for changes in in- 
terlayer spacing in pure fee and bcc crystals. In all cases, the semiempirical, theoretical 
techniques used rely either on input data (generally experimentally determined) or on cer- 
tain approximations for some of the variables of relevance. Necessarily, results will depend 
on such choices. Multilayer relaxations involve at best very small changes in position, and 
correspondingly, comparable changes in surface energy, whose minimization is the criterion 
used to determine the final interlayer spacings. Thus, the search for a minimum of the sur- 
face energy, as accurate as the minimization technique might be, will be strongly influenced 
by the two factors indicated above: the approximations used and the shallowness of the 
minimum in the surface energy surface resulting from small changes in the input parame- 
ters. As a consequence, to quote just one value for each of the changes in interlayer spacings 
as is ordinarily done, might not reflect the ambiguities in these calculations. In this paper 
we adopt a. different path: to each theoretical prediction, we will attach an estimate of the 
possible errors due to any of the reasons mentioned above. Although there is no certain 
way to determine such errors (after all, the predictions are, within their own framework, 
exact), we will see that changes on the order of 1 % in the surface energy can generate quite 
interesting variations in the relaxation schemes predicted. In particular, within the frame- 
work of ECT, such small changes in the surface energy can be easily obtained by changing 
any of the input parameters (lattice constant, cohesive energy, bulk modulus) by a similar 
amount, well below the usual experimental errors in the determination of such quantities. 
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To illustrate this issue, we will focus our attention on the surface structure of some 
fee pure metals (Al, Au, Cu and Ni). As can be seen in Tables 2*11 of ref. 1, previous 
theoretical and experimental studies show a wide spread in the predictions of the changes 
in interlayer spacings for the (100) and (110) surfaces. Even results obtained within the 
same theoretical technique (embedded atom method (EAM), EOT) do not agree with each 
other (due to different fitting procedures of the embedding function in the case of EAM 
and different input data in both cases). Although there is general qualitative agreement, 
regarding the contraction or expansion pattern found for successive layers, in some cases the 
absolute theoretical values show poor agreement with experimental results (see, for example, 
Al (100)). The ECT results (from refs. 1 and 13) also highlight this inconsistency. The 
difference between the values obtained in this work and those from previous applications of 
ECT is easily traceable to slightly different values of some of the input parameters. 

As mentioned above, in order to account for these and other ambiguities in the calcu- 
lation, we investigated the change in predicted relaxations due to small changes in the rigid 
surface energy. We thus defined ‘error bars’ in such way that all the intermediate values so 
obtained predict variations in surface energies within that tolerance. Needless to say, this 
range of values does not include all the possible sets (Adi 2 , Ad 23 ) that correspond to surface 
energies within the allowed values. It is interesting to note, however, that in most cases, all 
the experimental as well as theoretical predictions fall within the range of uncertainties in 
such procedure. 

It should be noted that when comparing our theoretical predictions with available 
experimental results, the error bars quoted in each case are similar in that the optimum 
relaxations are determined by minimization of some property by varying the input parame- 
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ters. To illustrate this point, we first discuss the surface energies and multilayer relaxations 
of the unreconstructed low-index surfaces of pure Al, Ni, Cu and Au crystals. In Table 
2 we display the ECT predictions for the surface energies and compare the results with 
typical experimental values for polycrystalline samples [15,16]. The agreement is excellent 
in all cases. We note that experimental values for the surface energies are for polycrystalline 
surfaces, thus could be strongly dominated by the predominant surface plane. 

In table 3 we compare results for the multilayer relaxations of the first two interlayer 
spacings for those cases for which recent experimental data is available [16-24]. Once again 
the agreement is excellent, as it was shown in previous ECT studies of surface structure 
[1]. The inclusion of the theoretical ‘error bar’, as mentioned above, allows for a better 
comparison with experiment as it shows that for most cases, small changes in the input 
parameters of the method suffice to account for the whole range of possible experimental 
results. The exceptions are Al(100) and Al(lll), where the outward relaxation of the surface 
layer has been attributed to an electron promotion effect [17]. Semiempirical methods 
(ECT, EAM. etc.), unless specifically designed to do so, do not generally allow for such 
fine electronic structure effects, thus it is not surprising that our results for Adi 2 in these 
cases predict surface layer contractions, even when the ‘error bar’ is taken into account. 
For completeness we also include results for the surface relaxation when only the top plane 
is allowed to relax, in order to single out correlations with subsequent interlayer spacing 
changes on the surface plane. Again, the agreement with available experimental data is 
very good in all cases. 

4. Multilayer relaxation of high-index surfaces 

We now discuss the application of ECT to the study of the surface structure of high- 
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index faces of fee (Al, Cu, Ni) and bee (Fe) metals. For each case, we computed the changes 
in interlayer spacing for the top six layers as well as the changes in registry. We follow the 
notation used in previous work on similar systems [2-7]. As discussed in Section 1 and 3, 
in this work we focus our attention not only on the absolute values of the relaxations, as 
obtained from experiment and predicted by the theory, but in the associated uncertainties 
as well. Although this last issue could be of importance when comparing the quality of the 
predictions, the main reason why it is highlighted here is to indicate the influence of external 
variables on the final results. We believe that because of the nature of the procedure used 
to perform the LEED analysis on the one hand, and the inherent uncertainties brought on 
by the simplifications adopted in designing semiempirical techniques, a thorough discussion 
of the results would not be complete if this issue was not appropriately addressed. In 
this spirit, figs. 1-4 summarize the ECT results for several fee and bcc systems and the 
corresponding LEED results for those cases for which experimental results are available. We 
have chosen this format for presenting the data for ease of comparison of experiment with 
theory. Fig. 1 displays results for Al faces. Al (210) and (331) display a similar behavior: 
large relaxation of the first interlayer spacing, followed by a comparable contraction of the 
second layer and an expansion of the third. The corresponding changes in registry are 
at best too small to definitely predict a trend in either case. As discussed in a previous 
application of ECT to Al (210) [12], the ECT results give a final configuration of higher 
symmetry and optimized ‘coverage’ of the space between atoms in the top layer. Whereas 
both theory and experiment agree in the case of registry changes in that the trends are the 
same and there is substantial overlap of the error bars, the agreement is less noticeable for 
the perpendicular relaxations, although the theoretical error bars clearly indicate that only 
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small changes of the input parameters are needed to improve the quantitative agreement. 
Fig. 2 shows similar results for fee (311) surfaces for Al, Cu and Ni. The trends found 
in the LEED results for Adi 2 are reproduced by ECT in all three cases. Although there 
are seemingly poorer results for the other relaxations, the principal point in this paper is 
that small shifts in the input parameters in the semiempirical method can bring trends 
into agreement. In the case of Ni(311), for which experimental values have been reported 
[7], there is good agreement with the associated registry changes. For completeness, we 
include the ECT results for Ni and Cu (210) (fig. 3), although there is no experimental 
data available for comparison. In spite of the differences in electronic structure, both metals 
display an almost identical behavior regarding the structure of the surface, which, together 
with the Al (210) results displayed in Fig. l.a, indicate a defined relaxation pattern for 

such fee faces. 

We conclude the presentation of results with a bcc system, Fe, for which several studies 
have been carried out. Fig. 4 displays LEED [3] and ECT results for Fe (210) and (310) 
surfaces. As expected, the percentage change in interplanar spacings is much larger for 
bcc metals than for fee, a feature clearly reproduced by ECT. The contraction-expansion 
pattern is also generally reproduced for both parallel and perpendicular relaxations. For 
completeness, we include the numerical results for all the cases studied in this paper in 
tables 4-7. Table 4 displays the results for the (210), (311) and (331) surfaces of Al where 
the error bars are related to changes of just 1 % in the surface energy. Table 5 and 6 show 
the results for the (311) and (210) faces of Ni and Cu, respectively. 

5. Conclusions 

In conclusion, in this paper we have used equivalent crystal theory to examine relax- 
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ations in higher index planes of Al, Cu, Ni and Fe where both perpendicular and parallel 
surface relaxations can occur and there is experimental data available for comparison. In 
addition, we attempt to raise questions regarding the nature of agreement between theory 
and experiment in that the success of theoretical agreement has been based on trends and 
absolute values of specific relaxations. Since semiempirical theoretical methods involve the 
use of experimental input parameters and conclusions are based on small changes in lattice 
geometry, we test the sensitivity of our predictions to small variations in input parameters. 
We find that in most cases, agreement involving trends can be greatly improved by small 
changes in these parameters and that careful consideration must be given to the method of 
making comparisons. 
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TABLE 1: COMPUTED CONSTANT AND EXPERIMENTAL INPUT FOR ECT 


[The constant p is 2n — 2, where n is the atomic principal quantum number, 1 (in A) is a 
scaling length and A (in A) is a screening parameter (see text). The constants A 3 and 
A 4 are dimensionless. A E (in eV) is the cohesive energy and a e (in A) the equilibrium 
lattice constant.] 


Element 

P 

l 

a 

A 

10 ~' 2 a 2 /d 

IO-'Aa/D 

10~ 4 D 

AE 

a e 

Al 

4 

0.336 

2.105 

0.944 

7.822 

2.104 

591.4 

3.34 

4.05 

Cu 

6 

0.272 

2.935 

0.765 

5.784 

2.530 

99.74 

3.50 

3.615 

Ni 

6 

0.270 

3.015 

0.759 

7.382 

2.793 

100.1 

4.435 

3.524 

Fe 

6 

0.277 

3.124 

0.770 

9.183 

1.887 

60.62 

4.29 

2.86 


TABLE 2: EXPERIMENTAL (Exp.) AND 
RELAXED ECT SURFACE ENERGIES 
OF Al, Cu, Ni AND Au 


[In ergs/cm 2 .] 


Technique 

Al 

Cu 

Ni 

Au 

Exp. [15] 



2270 

wm 

Exp. [16] 


Bran 

2380 

a 

Exp. [16] 

Brail 

llllJ 

2240 

1540 | 

ECT(IOO) 

1203 

2309 

2982 

wm 

ECT(llO) 

1284 

2373 

3073 

wm 

ECT(lll) 

856 

1767 

2274 

1136 | 


TABLE 3: SURFACE RELAXATIONS OF Al, Cu AND Ni AS PERCENTAGES OF THE 

BULK INTERPLANAR SPACINGS 


[The ECT A d 12 column displays results for relaxations of the top layer only while the ECT 
(two layers) columns display results for the case when the top two layers are allowed to 
relax.] 


Element 

Face 

Experiment 


ECT 

ECT 

(two- layers) 



Ad 12 

Ac?23 

Ref. 

Ad 12 

Ad 12 

Ac?23 


(100) 

+ 1.8 


[18] 

-4.68+1.62 

-5.05+1.58 

+3.35+0.80 

Al 

(110) 

-8.5+1.0 

+5.5+1. 1 

[19] 

-8.29+2.35 

-9.53+3.58 

+1.90+2.24 


(111) 

+1. 7+0.3 

+0.5+0. 7 

[20] 

-3.67+1.21 

-3.94+1.19 

+2.75+0.61 


(100) 

-3.2+0. 5 


[21] 

-3.53+1.68 

-3.82+1.68 

+2.48+0.85 

Ni 

(110) 

-9.0+1.0 

+3.5+1. 5 

[22] 

-6.32+2.44 

-6.55+3.63 

+0.34+2.24 


(111) 

-1.2+1.2 


[23] 

-2.89+1.29 

-3.10+1.25 

+2.12+0.63 



-2.1 

+0.45 

[24] 

-3.52+1.74 

-3.81+1.70 

+2.47+0.86 

Cu 

IBB 

-7.5+1.5 

+2.5+1.5 

[25] 

-6.31+2.46 

-6.51+3.83 

+0.29+2.44 


RSI 

-0.7+0. 5 


[26] 

-2.88+1.30 

-3.10+1.25 

+2.12+0.63 
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TABLE 4: PERPENDICULAR (Ad) AND PARAL- 
LEL (Aa) RELAXATIONS OF A1 (210), (311) 
AND (331) SURFACES EXPRESSED AS 
PERCENTAGES OF THE CORRE- 


SPONDING BULK SPACINGS 


Adi 2 

Ac?23 

Ari 34 

Ac /45 

Arise 

A1 (210) 

A1 (311) 

A1 (331) 

' -$.08±4.44 
-7.07±3.82 
+2.90±4.23 
-3.36±5.83 
+4.21±6.84 

- 8 . 88 ± 2.88 

-0.04±2.77 

-2.42±4.36 

+5.05±5.25 

-1.96±4.92 

-4.17±4.17 

-4.52±3.47 

+6.08±3.59 

-3.56±4.97 

+3.36±5.95 

Aai 2 

A(l23 

ACI 34 

AG 45 

A«56 

-0.20±2.42 

+0.02±2.60 

+0.79±2.90 

+0.04±3.49 

-0.47±4.42 

+0.87±2.67 

+0.53±3.19 

-0.81±4.09 

+2.27±3.86 

-0.37±4.96 

-2.65±2.26 

+0.02±2.63 

-0.08±3.08 

+1.87±3.64 

-1.41±3.76 


TABLE 5: PERPENDICULAR (Ad) AND 
PARALLEL (Aa) RELAXATIONS OF Ni 
(210) AND (811) SURFACES EX- 
PRESSED AS PERCENTAGES 
OF THE CORRESPONDING 
BULK SPACINGS 



Ni (210) 

Ni (311) 

Arii 2 

Ari23 

Ari 34 

Ari45 

Arise 

-4.54±5.04 

-4.96±4.15 

+1.06±4.37 

-1.96±5.76 

+3.17±6.73 

-5.57±3.26 
-0.77±3.00 
-1.61 ±4. 32 
+4.11±5.01 
-1.41±5.00 

A< 2 i 2 

A 023 
A 034 

Aa4s 

Aa56 

0.00±2.42 

+0.03±2.55 

+0.46±2.78 

+0.15±3.38 

-0.45±4.12 

+0.53±2.80 

+0.46±3.10 

-0.52±3.92 

+0.53±4.89 

-0.31±4.58 
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TABLE 6: PERPENDICULAR (Ad) AND 
PARALLEL (Aa) RELAXATIONS OF 
Cu (210) AND (311) SURFACES 
EXPRESSED AS PERCENT- 
AGES OF THE CORRE- 
SPONDING BULK 


SPACINGS 



Cu (210) 

Cu (311) 

Adu 

A^23 

Ac?34 

A^45 

A^56 

-4.48±5.07 

-4.91±4.17 

+0.96±4.38 

-2.03±5.80 

+3.20±6.80 

-5.54±3.27 

-0.81±2.99 

-1.75±4.33 

+4.32±5.05 

-1.60±5.02 

Aaj2 

Aa23 

Aa^ 

Aa 4 5 

Aa 56 

0.00±2.43 

+0.05±2.55 

+0.48±2.82 

+0.11±3.49 

-0.50±4.31 

+0.55±2.79 

+0.50±3.'12 

-0.54±4.00 

+0.50±5.08 

-0.30±4.73 


TABLE 7: PERPENDICULAR (Ad) AND 
PARALLEL (Aa) RELAXATIONS OF Fe 
(210) AND (310) SURFACES EX- 
PRESSED AS PERCENTAGES 
OF THE CORRESPONDING 


BULK SPACINGS 



Fe (210) 

Fe (310) 

Adi2 

Ad23 

Ad34 

Ad45 

Ad 5 6 

-7.65±7.06 

-2.81±5.48 

-4.53±4.84 

+0.06±5.14 

+1.00±4.30 

-30.45±1.99 

+15.5±1.83 

-19.82±1.91 

+11.49±1.81 

Aaj2 

A 023 

Aa34 

a 045 
A 056 

+2.37±2.12 

+3.10±2.05 

+1.47±2.70 

-0.75±2.70 

+0.40±3.23 

-2.36±1.47 

+6.36±1.67 

-4.82±1.92 

+2.60±1.93 
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Figure 3.— Theoretical values for the perpendicular (Ad) and 
parallel (Aa) relaxations for the (21 0) face of Cu (solid line) and 
Ni (dotted line), expressed as percentages of the corresponding 
bulk spaclngs. 


Fe ( 21 °) Fe (310) 



Figure 4.— Theoretical (solid line) and experimental (dotted line) [3] values for the perpendicular (Ad) and parallel (Aa) relaxations for the 
(a) (210) and (b) (310) faces of Fe, expressed as percentages of the corresponding bulk spaclngs. 
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